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Abstract 

We have made the first observation of B'^ ^Dg* Dg decays using 23.6 fb~^ of data recorded 
by the Belle experiment running on the T{5S) resonance. The branching fractions are measured 
to be B{B',^DtD~) = (1.03 tg;! tO;|)%, B{B',^Dt^Df) = (2.75 tg-.^^? ±0.69)%, and B(i?o^ 
D*^Dr) = (3.08tl;22 10;|)%; the sum is B{B^, ^ Di*^+ oi*^-) = {6.85t\fotlil)%. Assuming 
Bg ^ Ds* Ds* saturates decays to CP-even final states, the branching fraction determines the 
ratio ATg/ cos f, where AF^ is the difference in widths between the two B^-Bg mass eigenstates, 
and (/9 is a CP-violating weak phase. Taking CP violation to be negligibly small, we obtain 
Ar^/r^ = 0.147+0030 (stat.) t^^i (syst.), where T^ is the mean decay width. 



Decays of B^ mesons help elucidate the weak Cabibbo-Kobayashi-Maskawa structure of 
the Standard Model (SM). Because they are not produced in T(45') decays, B^ mesons are 
much less studied than their B^ and B^ counterparts. Most B^ data comes from the hadron 
collider experiments CDF and D0. Recently, another method to study B^ decays has been 
exploited: that of running an e^e^ collider at a center-of-mass (CM) energy corresponding 
to the T(5S') resonance, which subsequently decays to B^/^'B^^^ pairs. Both the CLEO [1| 
and Belle [21-^] Collaborations have used this method to measure inclusive and exclusive -B° 
decays. In this paper we use Lj^^^ = 23.6 fb~^ of data recorded by Belle at the T(5S') (i/i = 



10.87 GeV) to make the first observation (>5a) of 5°^ /}(*)+/}(*)- 
of the £)(*)+!)(*)- final state were made by the ALEPH ^ and D0 



5|. First measurements 
81] Collaborations using 



inclusive (fxpX and DfD~X samples. CDF [7| measured the single decay B'^ ^^ D'^ Dj . Here 
we exclusively reconstruct all three final states: DfD^, D*^Df, and D*^D*^ . 

These final states are expected to be predominantly CP-even J9|, and the (Cabibbo- 
favored) partial widths should dominate the difference in decay widths AF^'^ between the 
two Bg-Bg CP eigenstates J9|. This parameter equals AT^/ coscp, where AF^ is the width 
difference between the mass eigenstates, and </? is a CP- violating weak phase 13|. Thus 
the branching fraction gives a constraint in AF,-(/3 parameter space. Both of these param- 



10 



11| . The values favored by current 



eters can receive contributions from new physics 
measurements 12 1 differ somewhat from the SM prediction 10 1 . 

n n 

The Belle detector [IJ] running at the KEKB collider [15| includes a silicon vertex detec- 
tor, a central drift chamber, an array of aerogel threshold Cherenkov counters, time-of-flight 
scintillation counters, and an electromagnetic calorimeter. At the T(55') resonance, the 



e~^e — !■ bb cross section is a 



bb 



0.302 ± 0.014 nb 



na, 



and the fraction of T{5S) decays 



producing B^ mesons is f^ = 0.193 ±0.029 [16|. Three production modes are kinematically 
allowed: B^B^, B^B* or B*B^, and B*B*. In this analysis we use only the last (dominant) 



mode, for which the fraction is f^,^, = 0.901 Igolo 



n 



(|. The B* decays via B* -^ P,7, 
and the 7 is not reconstructed. Thus the number of B^B^ pairs used in this analysis is 

^B^B^ = Ant ■ ^bi ■ fs ■ fstB* = (1-24 ± 0.20) X 10^. 

We select 5° -^ D*+D*-, DfDf, and DfD- decays in which D+ -^ 07r+, i^^/s:+, 
K*^K~^, 0^"*", K^K*~^, and K*^K*'^. We require that charged tracks originate from near 
the e~^e~ interaction point. Charged kaons are selected by requiring that a kaon likelihood 
variable based on dE/dx measured in the central drift chamber and information from the 



aerogel threshold Cherenkov counters and time-of-flight scintillation counters be >0.60; this 
requirement is ~90% efficient and has a vr^ misidentification rate of ~10%. Tracks having 
kaon likelihood < 0.60 are identified as vr^. Neutral K'g candidates are reconstructed from 



vr ' TX 



pairs having an invariant mass within 10 MeV/c^ of the K^ mass 16| and satisfying 
loose requirements on the decay vertex position [l7| . The momentum of tracks (except the 
vr^ from K^ decay) must be > 100 MeV/c. 

Neutral tt^ candidates are reconstructed from 77 pairs having an invariant mass within 
15 MeV/c^ of the vr" mass. The photons must have a laboratory energy greater than 
100 MeV. Neutral K*'^ (charged K*^) candidates are reconstructed from a K~ (if^) and 
TT"*" having an invariant mass within 50 MeV/c^ of M^.o (M^,+ ). Neutral (f) (charged p^) 
candidates are reconstructed from a K^K~ (tt+vt") pair having an invariant mass within 
12 MeV/c2 (100 MeV/c2) of M^ (M^+). 

The invariant mass windows used for D^ candidates are 10 MeV/c^ (2.5 — 3.2cr) for the 
three final states containing K* candidates, 20 MeV/c^ (l-^o") for 0p^, and 15 MeV/c^ (^ 
4. Oct) for the remaining two modes. For the three vector-pseudoscalar final states, we impose 
a loose requirement on the helicity angle ^^ob which is the angle between the momentum of 
the charged daughter of the vector particle and the direction opposite the Dg momentum 
in the rest frame of the vector particle. We require | cos^j^^jjl > 0.20, which retains 99% of 
signal decays and rejects 18% of remaining background. 

To reconstruct D*^ — )■ -D^7 decays, we pair D^ candidates with photon candidates and 
require that the mass difference M^+ — M^+ be within 12.0 MeV/c^ of the nominal value 
(143.8 MeV/c^), where D^ denotes the reconstructed -D+ candidate. This requirement (and 



that for TZ discussed below) is determined by optimizing a figure-of- merit S/y/S + B, where 
5* is the expected signal based on Monte Carlo (MC) simulation and B is the expected 
background as estimated from a data sideband. We require that the photon energy in the 
CM system be greater than 50 MeV, and that the energy deposited in the central 3x3 array 
of cells of the electromagnetic calorimeter cluster be at least 85% of the energy deposited in 
the central 5x5 array of cells. 

Signal Bg decays are reconstructed from D^*^D^g*^ pairs using two quantities: the beam- 



energy-constrained mass My^^ = J E^^^^ — p\, and the energy difference Ai5 = E^ — -Ebcam; 
where p^ and E^ are the reconstructed momentum and energy of the 5°, and -Ebcam is the 
beam energy. These quantities are evaluated in the e"'"e~ CM frame. When the i?° is not 



fully reconstructed, e.g., due to losing the 7 from D*+ — )• -D+7, AE is shifted lower but 
M-^^ remains almost unchanged. We determine our signal yields by fitting events in the 
region 5.20 GeV/c^ < M^, < 5.45 GeV/c^ and -0.15 GeV < AE < 0.10 GeV. The modes 
T(5S') —7- i?^i?g, B^B* and B*B* are well-separated in My^^-AE space. We see no evidence 
for B^Bg and B^B* and thus do not fit for them. The expected yields based on Ref. |j] are 
less than one event for each of DfD~, D*^Df, D*'^D*~ final states. 

Approximately half of the events have multiple fi^ — j-D^*^"*"/}^*^" candidates, which usu- 
ally arise from low momentum 7's produced from 7r° decays. For these events we select the 
candidate that minimizes the quantity 

1 r 



X' = TTT^l^JE [{MD.-M^J/a,,]' + Y. [{AM - AM)/a^,,f] , (1) 



,2 _ 

(2 + AT) 



where AM = M^* — M^ ; M^ and AM are reconstructed quantities; a^,,^ and a^^j are the 
uncertainties on M^ and AM; and the summations run over the two Df daughters and 
possible D*~^ daughters (A^ = 0, 1,2) of a 5° candidate. According to the MC simulation, 
this criterion selects the correct B^ candidate 85%, 76%, and 75% of the time, respectively, 
for DfD-, D*^Df, and D*+D*- final states. 

We reject background from e~^e~ -^ qq {q = u, d, s, c) continuum events based on event 
topology: qq events tend to be colhmated, while B^^Bi^^ events tend to be spherical. We 
distinguish these topologies using a Fisher discriminant based on a set of modified Fox- 



Wolfram moments 18|]. This discriminant is used to calculate a likelihood Cs {Cqq) for an 
event assuming the event is signal {qq background). We form the ratio IZ = Cs/{Cs + Cq^ 
and require TZ > 0.20. This selection is 95% efficient for signal decays and removes > 80% 
of qq background. 

The remaining background consists of T {5S) —)■ B^/W^*^ —)■ D'^ X , T{5S)^BBX (where 
bb hadronizes into 5°, 5°, or B^), and B^ -^ Dfj{2317)Di*\ B^ -^ Dfj{24:60)Di*\ and 
Bg — )■ DfDjTT^ decays. The last three processes peak at negative values of AE, and their 
yields are estimated using analogous B^ — )■ D^jD^*^ branching fractions. The total yields 
for all backgrounds within an M^^^-AE signal region spanning 3a in (M^^^, AE) resolution 
are 0.25 ± 0.03, 0.25 ± 0.06, and 0.15 ± 0.13 events, respectively, for B^^DfD;, D*^Df, 
and D*^D*~~ decays. To check our background estimates, we count events in the sideband 
region M^^ < 5.375 GeV/c^ and find reasonable agreement with the yields predicted from 
MC simulation. All selection criteria are finalized before looking at events in the signal 
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regions. The final event samples are shown in Fig. [H 

To measure the signal yields, we perform a two-dimensional extended unbinned 
maximum-likelihood fit to the M-^^^-^E distributions. For each sample, we include prob- 
ability density functions (PDFs) for signal and qq, 5(*)fiW ^ £)+x, and T(55) -^ BBX 
backgrounds. As these backgrounds have similar Mj^^, AE' shapes, we use a single PDF for 



them, taken to be an ARGUS function [20| for M^^ and a third-order Chebyshev polynomial 
for ^E. All shape parameters are taken from MC simulation. Other backgrounds are very 
small and considered only when evaluating systematic uncertainties. 

The signal PDFs have three components: correctly reconstructed (CR) decays; "wrong 
combination" (WC) decays in which a non-signal track or photon is included in place of a 
true daughter track or photon; and "cross-feed" (CF) decays in which a D*^Df or D*'^D*~ is 
reconstructed as DfDj or D*g^DJ , respectively, or else a DfDj or D*g^DJ is reconstructed 
as D*'^Df or Z}*+D*~. In the former (latter) case the signal decay has lost (gained) a 
photon, and AE is typically shifted lower (higher) by 100-150 MeV. The PDF for CR 
events is modeled with a single Gaussian for M^^ and a double-Gaussian with common 
mean for AE. The means and widths are taken from MC simulation and calibrated using 
5°^D(*)-7r+ and B^^D^*^+D- control samples. The PDFs for WC and CF events are 



modeled from MC simulation using non-parametric PDFs with Kernel Estimation [19|. The 
fractions of WC and CF-down events are also taken from MC simulation. The fractions 
of CF-up events are difficult to simulate and thus floated in the fit, as the extraneous 7 
usually originates from a B decay chain and many B, B^ partial widths are unmeasured. As 
the CF-down fractions are fixed, the three distributions {DfDj, D'^'^Df, and D*+Z}*") are 



fitted simultaneously J2l|. The CF fractions are typically 0.1-0.4. 

The fit results are listed in Table HI and projections of the fit are shown in Fig. |2j The 
branching fraction for channel i is calculated as B^ = Yj{e\jQ-N^ -^ -2), where X is the fitted 



CR yield, and e\^(j is the MC efficiency with intermediate branching fractions 16| included. 



The efficiencies £5\//c include small correction factors to account for differences between MC 
simulation and data for kaon identification. Inserting all values gives the branching fractions 



listed in Table [H The statistical significance is calculated as J —2\\i{Cq/ C^^x), where £0 
and £max are the values of the likelihood function when the signal yield Y^ is fixed to zero 
and when it is the fitted value, respectively. We include systematic uncertainty in the 
significance by smearing the likelihood function by a Gaussian having a width equal to the 



total systematic error related to the signal yield. 

TABLE I: Signal yields (Y), efficiencies including intermediate branching fractions (e), branching 
fractions (B), and signal significance (S) including systematic uncertainty. The first error listed 
is statistical, the second is from systematics due to the analysis procedure, and the third is from 
systematics due to external inputs. 



Mode 


Y 

(events) 


e 
(xlO-4) 




B 

(%) 




S 


DtD- 


o c +3.2 

o-t'-2.6 


3.31 


1.03 


+0.39 +0.15 
-0.32 -0.13 


±0.21 


6.2 


Dl^Df 


9.2 +ij 


1.35 


2.75 1|]:?? ±0.40 ±0.56 


6.6 


Dl+Dl- 


4-911:? 


0.643 


3.08 


+1.22 +0.57 
-1.04-0.58 


±0.63 


3.1 


Sum 


22.6 If] 




6.85 


+1.53 1111 +1.40 
-1.30 =■= ^-^^ -1.41 





The systematic errors are listed in Table [TTl The error due to PDF shapes is evaluated by 
varying shape parameters by ±lcr for backgrounds, and by trying different parameterizations 
for the WC and CF components. The systematic error for the fixed CF-down fractions is 
evaluated by fitting a i?^ — )■ D~tt^ control sample and comparing the fraction of 5° — )■ D*~n^ 
contamination with that predicted by MC simulation. The difference is taken as the range 
of variation for the CF-down fractions. We vary the fractions over this range and take 
the r.m.s. variation in B as the systematic error. The error due to the fixed WC fractions 
is evaluated in the same way, but the range of variation for the WC fractions is ±20%. 
The uncertainty due to K^ identification depends on momentum but is ~2.5% per track; 
as our final states typically have four charged kaons, this error is 10-11%. The error due 
to tracking efficiency is 1.0% per track. Uncertainty due to the unknown B^ — )■ D*'^D*~ 
longitudinal polarization fraction (/j^) affects all three modes due to the CF components. 
For our nominal result, we take fi^ to be the world average (WA) value for the analogous 
spectator decay B^-^D^^D*^: 0.52 ± 0.05 [16|. The systematic error is taken as the change 
in B when /^ is varied by twice the error on the WA value. Significant uncertainties arise 
from D^ branching fractions, O'-^i^s)^ fsi ^"^^ fn,-n*i which are external factors that should 
be measured more precisely in the future. We list separately the systematic error due to 
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FIG. 1: M^^ vs. AE scatter plots. The signal ellipses correspond to 3a in resolution for T(5S')— )• 
B*Bl decays; the number of candidates within the ellipses is listed. The top, middle, and bottom 
plots correspond to B'^^^DfDj, B°^D*^Df, and B^^D*+D*-, respectively. 
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FIG. 2: Mj^j, and Ai? projections of the fit result. The rows correspond to B^ — >■ DfDg (top), 
B'^^D*^Df (middle),and B'^^D*+D*~ (bottom). The red dashed (dotted) curves show RC 
(WC) signal, the green and blue dash-dotted curves show CF signal, the grey solid curve shows 
background, and the black solid curves show the total. 

these factors in Table [11 

In the limit m^^ — t- oo while (m^ — 2m^) — )■ 0, the b — )■ ccs process saturates the decay 



width [22|. If also the number of colors A^^ -)■ oo, then B^ -)■ /}*+/}*", -0*^^-0^ (along 
with DfD-) are CP even (+), and r[B%CP+) -^ Di*Wf^ saturates AFf^ ^. This 



gives the relationship 2-6(5° ^DW+^W-) = (AFf^/2)[(l + cos^)/Fj 



cos if) /Tjj], 
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TABLE II: Systematic errors (%). The first twelve sources affect the signal yield and thus the 
signal significance. 



Source 


Dt 


Ds 


Dt^Df 


Dt+nr 




+a 


—a 


+a —a 


+a 


—a 


CR PDF Shape 


0.8 


0.8 


0.3 0.3 


0.5 


0.4 


Background PDF 


1.1 


1.3 


1.9 2.0 


3.0 


6.1 


WC+CF PDF 


0.3 


0.3 


1.5 1.5 


4.4 


4.4 


WC/CF Fractions 


0.2 


0.2 


5.0 5.0 


8.7 


8.7 


TZ Requirement [qq suppr.) 


1.8 


1.8 


1.8 1.8 


1.8 


1.8 


Best Candidate Selection 


6.9 


0.0 


2.2 0.0 


2.2 


0.0 


K^ Identification 


10.1 


10.1 


10.6 10.6 


10.9 


10.9 


K^ Reconstruction 


2.1 


2.1 


2.1 2.1 


2.2 


2.2 


TT^ Reconstruction 


1.1 


1.1 


1.1 1.1 


1.0 


1.0 


7 


- 


- 


3.8 3.8 


7.6 


7.6 


Tracking 


6.2 


6.2 


6.2 6.2 


6.2 


6.2 


Polarization 


0.2 


0.0 


0.8 0.5 


0.7 


0.3 


MC Statistics for e 


1.1 


1.1 


0.9 0.8 


1.0 


1.0 


Ds Branching Fractions 


12.4 


12.4 


12.4 12.4 


12.5 


12.5 


Luminosity 






±1.3 






^T{bS) 






±4.6 






fs 






±15 






^b*b: 






+4.2 
-4.4 






Total 


24.9 


24.0 


25.1 25.1 


27.5 


28.0 



where F^ j^ are the decay widths of the hght and heavy mass eigenstates [13|] . Substituting 



LM 



Q. 



r±Ar^/2 and AT^ = AT ^/ cos if [13| allows one to use the branching fraction B to 



constrain AF^ and (f. If CP violation is negligible, then cosc/?~l and the above expression 
can be inverted to give AF^/F^ = 2B/{1 — B). Inserting B from Table [T] yields 

AF„ 



0.147 



+0.036 +0.042 
-0.030 -0.041 ! 



(2) 



where the first error is statistical and the second is systematic. This result is 1.3o" higher than 
that of Ref. (8| but consistent with the theory prediction 10|. There is theoretical uncertainty 
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arising from the CP-odd component in B^ ^Dg'^Df, D*'^D*~ and contributions from other 
two-body final states; the effect upon AF^/F^ is estimated in Ref. 9| to be ±3%. This is 
much smaller than the statistical/systematic errors on our measurement, but there rnay be 
additional contributions coming from three-body final states, which are neglected in 9|. 

In summary, we have measured the branching fractions for 5° — t- D^*^~^ D^J*^~ using 
e"'"e~ data taken at the T{5S) resonance. Our results constitute the first observation of 
B^^D*'^Df (6.6cr significance) and provide the first evidence for B^^D*^D*^ (3.1(T sig- 
nificance). We use these measurements to determine the B'^-B^ decay width difference AF^ 
with improved fractional precision. 
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